The association between Toxoplasma infection and mortality: the NHANES epidemiologic follow-up study

Background Toxoplasma gondii has been reported to be associated with higher mortality in patients with schizophrenia. This study aimed to explore the relationship between T. gondii infection and 25-year mortality based on data from the Third National Health and Nutrition Examination Survey (NHANES III) database. Methods Cases with serum T. gondii antibody test results were included in this study and the corresponding mortality dataset was obtained from the US National Center for Health Statistics (NCHS). Propensity score matching (PSM) was used to match age and sex between groups. The Cox proportional hazards model was used to evaluate the effect of T. gondii infection on mortality. Results A total of 14,181 cases were included in the analysis, of which 3831 (27.0%) were seropositive for T. gondii antibody. The median follow-up time of the whole cohort was 22.5 (interquartile range 16.3, 24.5) years. A total of 5082 deaths were observed in this cohort, a mortality rate of 35.8%. All-cause mortality was significantly higher in the seropositive group than in the seronegative group (50.0% vs 30.6%, P < 0.001). Kaplan–Meier analysis showed a significant difference in the survival time between two groups before and after PSM. Multivariate analysis showed that T. gondii infection was independently associated with higher all-cause mortality after adjusting for potential confounders. Conclusions Toxoplasma gondii infection is associated with higher mortality in general population. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s13071-022-05398-1.


Background
Zoonoses continue to be a major public health concern [1]. Toxoplasma gondii is a common protozoan parasite that infects approximately 30% of the world's population [2,3]. In healthy individuals, T. gondii infection usually results in mild illness most of the time, but it has the potential to cause severe consequences. Toxoplasmosis is a common opportunistic infection in immunocompromised patients and results in high mortality in this patient population [2,4,5].
Toxoplasmosis has been linked to a number of diseases, including neuropsychiatric diseases, such as primary neuropathies, behavioral and psychiatric disorders [6,7], as well as some related metabolic diseases, such as non-alcoholic fatty liver disease [8] and diabetes [9,10]. Toxoplasma gondii infection is also associated with a worse prognosis of all-cause mortality in individuals with these diseases; for example, it has been shown to be associated with higher all-cause mortality rate in patients with schizophrenia [11,12]. However, it remains unclear whether T. gondii infection actually affects the all-cause mortality rate in the general population.
The present study aimed to explore the relationship between T. gondii infection and 25-year mortality in the general population based on data retrieved from a nationally representative longitudinal database.

Study population
The study population was from the Third National Health and Nutrition Examination Survey (NHANES III) 1988-1994, which was a periodic survey conducted by the US Centers for Disease Control and Prevention. The study data are freely available for use by researchers from all over the world. We selected cases with serum T. gondii antibody test results for the present study. All individuals were followed up until December 2015, and the data on survival status and death dates were collected. The corresponding mortality dataset was obtained from the US National Center for Health Statistics (NCHS), which is also freely accessible. All NHANES protocols had been approved by the NCHS ethics review board, and participants or their proxies had provided informed consent prior to participation. As the present study only utilized anonymized data available to the public, no additional ethics approval was required.

Definitions
The presence of T. gondii antibody was tested using the indirect enzyme immunoassay method. The diagnostic threshold for T. gondii infection was 7 IU/ml as described in the original report. Samples with results ≥ 7 IU/ml were considered to indicate T. gondii infection. The entire dataset and details on the methods used can be retrieved from the NHANES website (https:// www. cdc. gov/ nchs/ nhanes/ about_ nhanes. htm).
Hypertension was defined as a systolic blood pressure ≥ 135 mmHg or diastolic blood pressure ≥ 85 mmHg. Patients with a significant history of hypertension and/or undertreatment of hypertension were also included.
People with type 2 diabetes were defined as having a history of diabetes, currently using insulin or oral hypoglycemic agents, with a fasting blood glucose (FBG) ≥ 7.0 mmol/l or glycated hemoglobin (HbA1c) ≥ 6.5% or 2-h postprandial glucose ≥ 11.0 mmol/l. Body mass index (BMI, kg/m 2 ) was defined as the weight in kilograms divided by the square of the height in meters.
Income was measured based on the poverty income ratio, an index reflecting the ratio of household income to the household poverty level determined by area of residence and household size. Low family income was defined as a ratio below 1.
In NHANES III, data on education was collected as a numerical variable (number of years of education), ranging from zero to 17 years. Low educational level was defined as an educational level of lower than the 12th grade.

Statistical analysis
Categorical variables were expressed as percentages. Continuous variables were expressed as means ± standard deviation. The Student t-test (for variables normally distributed), the Mann-Whitney U-test (for variables non-normally distributed) and the Chi-square test (for categorical variables) were used to compare the differences between the seropositive and seronegative groups. We used propensity score matching (PSM) to select ageand sex-matched participants between the two groups. Patients were matched 1:1 based on their propensity scores, and the match tolerance value was 0.0001. The Cox proportional hazards model was used to evaluate the effect of T. gondii infection on all-cause mortality. All tests were two-tailed and a P-value < 0.05 was considered to be statistically significant. All analysis was conducted in R 3.6.2 (https:// www.r-proje ct. org/).

Baseline characteristics of the study population
A total of 19,598 participants with follow-up data were screened for this study (Fig. 1). After excluding participants without T. gondii antibody test (N = 4,085), 1332 without other key data, 14,181 cases were included in final analysis. Among them 6,672 (47.0%) were male and the average age was 46.3 ± 19.5 years old. A total of 2128 (15.0%) had diabetes and 5859 (41.3%) had hypertension. A total of 3831 (27.0%) cases were seropositive in T. gondii antibody. More details were shown in Table 1.

Comparison between the T. gondii-seropositive group and the T. gondii-seronegative group before and after PSM
According to the results of the T. gondii antibody test, participants were divided into a T. gondii-seropositive group and a T. gondii-seronegative group. The comparison of the baseline characteristics before PSM are shown in Table 1. Patients in the seropositive group were significantly older than their seronegative counterparts (mean ± SD, 53.7 ± 19.7 vs. 43.6 ± 18.7; P < 0.001). There was a higher proportion of males in the seropositive group than in the seronegative group (49.9% vs. 46.0%; P < 0.001). Diabetes, hypertension and metabolic derangements were more common in the seropositive group. Table 1 Comparison of Toxoplasma gondii-seropositive and -seronegative groups before propensity score matching Values in table are presented as the mean ± standard deviation (SD) or as the median with the interquartile range in square brackets, unless indicated otherwise ALT Alanine aminotransferase, AST aspartate aminotransferase, BMI body mass index, BUN blood urea nitrogen, CRP C-reactive protein, FBG fasting blood glucose, HbA1c glycosylated hemoglobin, BUN blood urea nitrogen.

Variables
Total T. gondii antibody P  As age and sex are factors closely associated with life expectancy, and a significant difference in these factors was found between the seropositive and seronegative group, we used PSM to match age and sex between two groups. A total of 3806 pairs of cases were selected in which age and sex were comparable after PSM (Table 2). It was noteworthy that the severity of metabolic derangements, the proportion of patients with diabetes and BMI level remained significantly different between two groups even after matching for age and sex.

All-cause mortality in T. gondii-seropositive and -seronegative groups
The median follow-up time of the whole cohort was 22.5 (interquartile range 16.3, 24.5) years. A total of 5082 deaths were observed in this cohort, which is a mortality rate of 35.8%. All-cause mortality was significantly higher in the seropositive group than in the seronegative group (50.0% vs. 30.6%; P < 0.001). Kaplan-Meier analysis showed a significant difference in the survival time between two groups (Fig. 2a). After PSM, the difference in the mortality rate between the seropositive group and seronegative group was attenuated but remained statistically significant (49.7% vs. 47.3; P = 0.037). The same result was obtained in the Kaplan-Meier analysis (P = 0.014; Fig. 2b).

Cox multivariate analysis for the impact of T. gondii infection on all-cause mortality after PSM
We used Cox analysis to explore the correlation between T. gondii infection and all-cause mortality in the PSM population. The results are shown in Table 3. In the univariate analysis, T. gondii antibody seropositivity, race, sex, age, low educational level, hypertension, diabetes, waist circumference, related metabolic parameters and liver function and renal function test results were significantly different between the survival group and mortality group. Variables with a P-value < 0.05 in Table 3 were considered to be candidates for inclusion in the multivariate regression. Collinearity analysis showed a strong collinearity between BMI and waist, FBG and HbA1c, AST and ALT, and creatinine and BUN. For the variables with collinearity, we used the receiver operating characteristic (ROC) analysis to select the variable with a higher area under the ROC curve. The results showed that waist, HbA1c, AST and BUN were the most suitable for inclusion in the multivariate Cox regression analysis.
In order to explore the association between T. gondii infection and mortality, we established three models with the combination of different variables to adjust for different confounding factors. Model 1 was adjusted for race, sex and   Fig. 2c).

The difference in risk factors for mortality in the T. gondii-seropositive and T. gondii-seronegative groups
The results of the Cox multivariate analysis in the T. gondii-seropositive and T. gondii-seronegative groups are shown in Additional file 1: Table S1. The effect of most factors was similar in the two groups, with the exception of CRP. CRP was significantly associated with the risk of mortality in the T. gondii-seropositive group (HR 1.070, 95% CI 1.020-1.121, P = 0.005), while no statistical significance was found in the seronegative group (P = 0.121).

Discussion
The present study used data from NHANES III to assess the longitudinal relationship between T. gondii infection at baseline and subsequent mortality after adjusting for relevant confounding factors. The main finding of this prospective study is that T. gondii antibody positivity was associated with increased all-cause mortality. Toxoplasmosis is generally regarded to be a threat to immunocompromised patients or those with congenital toxoplasmosis [2,13]. Several studies have also shown a higher risk for mortality from natural causes in persons with schizophrenia infected with T. gondii [4,11,12,14]. Our study is the first to show that T. gondii infection is associated with higher mortality in the general population.
Several possible hypotheses might explain this finding. First, T. gondii infections are associated with elevated biomarkers of chronic inflammation, such as CRP [15]. Systemic chronic inflammation can lead to several diseases and cause disability and mortality worldwide [16]. The results of the present study show that the CRP level was associated with higher mortality in the T. gondiiseropositive group, indicating the possible role of chronic inflammation in an increased mortality. Second, T. gondii can influence host metabolism of fatty acids, lipids and energy in the liver [17]. Higher seroprevalence of T. gondii antibody has been found in patients with non-alcoholic fatty liver disease [8], type 2 diabetes [9] and cardiovascular disease burden [3]. The higher risk of mortality in the toxoplasmosis population might be partially attributed to the increased metabolic morbidity following T. gondii infection. Third, toxoplasmosis has been showed to be linked to some neuropsychiatric disorders [18] and a higher risk of traffic accidents and suicide attempts [19,20]. The mechanism underlying the impact of toxoplasmosis on emotional and mental problems includes proinflammatory immune response to infection and endocrine and neurotransmitter dysregulation [19].
This study was based on a general population survey database and had a larger sample size and longer follow-up duration than other similar studies. Although the data are strong due to the number and quality of the data, as well as the multiple confounders we have tried to control, there are many possible explanations of this statistically significant association. Toxoplasma gondii infection has been shown to be associated with deregulated metabolic functions [21,22] and cardiovascular disease burden [3]. The factors that increase metabolic diseases can also increase the risk for toxoplasmosis acquisition. We have attempted to control various known confounding factors, including metabolic profiles, socioeconomic status and the medical history, but even with very good multivariate models, using this type of data it is impossible to separate which factor is really influencing the outcome. Notably, the effect of toxoplasmosis on mortality was not as strong as that of diabetes, hypertension and sex. Moreover, as shown in the Fig. 2, the entire effect of toxoplasmosis attenuated when age and sex were controlled. Therefore, we are unable to rule out the possibility that the rest of the effect was mediated by other, yet unknown confounding factor(s). A well-designed prospective study involving people with metabolic risk may possibly answer this question.

Conclusions
In conclusion, in this epidemiologic follow-up study, toxoplasmosis was found to be associated with higher all-cause mortality in the general population. The results suggest that toxoplasmosis remains a neglected parasitic infection which requires more public health action.